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Abstract 

We report many-body perturbation theory calculations of excited-state properties of distorted 
1-T diamond-chain monolayer rhenium disulfide (ReS 2 ) and diselenide (ReSe 2 ). Electronic self¬ 
energy substantially enhances their quasiparticle band gaps and, surprisingly, converts monolayer 
ReSe 2 to a direct-gap semiconductor, which was, however, regarded to be an indirect one by 
density-functional-theory calculations. Their optical absorption spectra are dictated by strongly 
bound excitons. Unlike hexagonal structures, the lowest-energy bright exciton of distorted 1-T 
ReS 2 exhibits a perfect figure-8 shape polarization dependence but those of ReSe 2 only exhibit a 
partial polarization dependence, which results from two nearly-degenerated bright excitons whose 
polarization preferences are not aligned. Our first-principles calculations are in excellent agreement 
with experiments and pave the way for optoelectronic applications. 
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I. INTRODUCTION 


To overcome the zero band gap of graphene^, two-dimensional (2D) transition-metal 
dichalcogenides (TMDCs) semicondnctors have attracted signihcant attention!^*-. Thanks 
to the very thin atomic thickness (a few A), pristine interfaces withont ont-of-plane dan¬ 
gling bonds, and considerable band gaps (1-2 , TMDCs are potentially benehcial 

for eliminating short channel effects^, lowering interface state densities, redncing snrface 
ronghness scattering!^, and low-power digital applications!^^—. However, becanse of inter¬ 
layer interactions and variations of screening, the electronic strnctnres of many hexagonal 
TMDCs nndergo snbstantial variations with different stacking layer nnmbers. For example, 
a crossover from direct band gaps in monolayers to indirect band gaps in mnltilayers has 
been observed in MoS 2 ^‘!‘!!, limiting the applicability of TMDCs in optoelectronic devices. 

Layered crystals of rheninm disnlhde (ReS 2 ) and diselenide (ReSe 2 ) are a new family of 
2D TMDCs semicondnctors, which have been snccessfnlly fabricated recently!^'—. Unlike 
hexagonal TMDCs, ReX 2 (X= S, Se) crystallizes in a distorted 1-T diamond-chain strnctnre 
with the triclinic symmetry, as a resnlt of charge deconpling from an extra valence electron 
of Re atoms!^^— . This strnctnral distortion leads to a mnch weaker interlayer conpling. 
Conseqnently, the band renormalization is absent and bnlk ReX 2 behaves as electronically 
and vibrationally deconpled monolayersl^i!^. Snch a vanishing interlayer conpling in ReX 2 
strnctnres enables probing 2D-like systems withont the need of monolayer or few layers, 
overcoming the challenge of preparing large-area, single-crystal monolayers!^. Fnrthermore, 
the strnctnre distortion of lT-ReX 2 makes it exhibit nniqne anisotropic optical, electrical, 
and mechanical properties!^!^, with considerable interests for varions applications in po¬ 
larization controller, liqnid crystal displays, 3D visnalization techniqnes, (bio)dermatology, 
and in optical qnantnm compnters^ii^^. 

Despite these nniqne properties and potential applications of distorted 1-T diamond- 
chain ReX 2 strnctnres, we have very limited knowledge abont their fnndamental excited- 
state properties, snch as qnasiparticle band gaps, optical spectra, and excitonic effects. In 
particnlar, currently available experimental measurements are diverse due to sample qualities 
and unavoidable environment effects. In this sense, a reliable and parameter-free calculation 
to accurately capture essential electron-electron (e-e) and electron-hole {e-h) interactions 
in excited-state properties is highly desirable for understanding fundamental physics and 
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motivating applications of few-layer ReX2. 

In this work, we employ the hrst-principles GW-Bethe-Salpeter Equation (BSE) approach 
to study two typical distorted 1 -T diamond-chain TMDCs, i.e., monolayer ReS2 and ReSe2. 
We observe signihcant many-electron effects that dictate their excited-state properties. Self¬ 
energy corrections not only enlarge band gaps but also convert suspended monolayer ReSe2 
to be a direct-gap semiconductor. Excitonic effects are prominent and the BSE-calculated 
absorption spectra are in excellent agreement with measurements. In particular, because 
of the reduced triclinic symmetry, monolayer ReS2 and ReSe2 possess anisotropic optical 
responses within the near-infrared frequency regime. These excited-state properties within 
the interesting frequency regime may be useful for optoelectronic applications. 

The remainder of this paper is organized as follows: In Sec. II, we introduce the atomic 
structures of monolayer ReS2 and ReSe2, and our computational approaches. In Sec. Ill, 
the quasiparticle band energies, band gaps, and the transition from indirect to a direct band 
gap, are presented. In Sec. IV, we discuss the optical absorption spectra and excitonic 
effects. In Sec. V, we focus on the anisotropic optical responses of excitons. In Sec. VI, we 
summarize our studies and conclusion. 


II. ATOMIC STRUCTURE AND COMPUTATIONAL METHODS 

The ball-stick structures of monolayer distorted 1 -T diamond-chain ReS2 and ReSe2 are 
presented in Fig. [T] The unit cell and unit vectors are marked by purple-colored vectors. 
Each Re atom is clustering of “diamond chains” composed of the zigzag shape along the a 
direction as shown by red-colored lines in Fig. [T](a)^‘^^— . The side view of monolayer ReS2 
or ReSe2 is presented in Fig. [ 11 (b). The special feature is that, unlike hexagonal TMDCs, 
the S or Se atoms of distorted 1 -T diamond-chain structures are not all in the same plane, 
substantially lowering the structure symmetry. The hrst Brillouin zone (BZ) is plotted in 
Fig. HKc), which is a hexagon but with unequal side lengths as a result of the distorted 
atomic structure. Consequently, the Kl, K 2 , and K 3 points are no longer equivalent. 

We fully relax the structures according to the force and stress calculated by DFT with 
the Perdew, Burke, and Ernzerhof (PBE) functional^S, using the Quantum Espresso 
package^!.. The ground-state wave functions and eigenvalues are obtained from the DFT/PBE 
with norm-conserving pseudopotentials^. The plane-wave basis is set with a cutoff energy 
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of 60 Ry with a 8 x 8 x 1 fc-point grid. A vacuum space between neighboring layers is set 
to be more than 25 A to avoid the interactions between layers. Based on these parameters, 
our calculated lattice constant a is 6.43 (6.78) A and the lattice constant b is 6.53 (6.66) A 
for ReS 2 and ReSe 2 , respectively. They are consistent with previous calculationsi^i^. 

The quasiparticle energies and band gaps are calculated by the GW approximation within 
the general plasmon pole (GPP) model^. The involved unoccupied conduction band num¬ 
ber for calculating the dielectric function and self-energy is about ten times of the occupied 
valence band number. In solving the BSE, we use a hner /c-point grid of 32 x 32 x 1 for 
converged excitonic states^. All the GW-BSE calculations are performed with the Berke- 
leyGW code^ including the slab coulomb truncation scheme to mimic suspended monolayer 
structures^^i^. For optical absorption spectra, only the incident light polarized parallel with 
the plane is considered due to the depolarization effect^*^. 


III. QUASIPARTICLE ENERGIES AND BAND GAP 

The DFT calculated band structures of monolayer ReS 2 and ReSe 2 are presented in 
Fig. [21 That with spin-orbit coupling (SOG) included is presented in the appendix. As 
shown in Fig. [2](a), monolayer ReS 2 exhibits a 1.36 eV direct band gap at the P point. 
This is consistent with a recent work— reporting a value of 1.43 eV. On the other hand, 
at the DFT level, monolayer ReSe 2 is observed to be an indirect-gap semiconductor; the 
conduction band minimum (GBM) is located at the P point while the valence band maximum 
(VBM) is slightly away from the P point, as shown in Fig. [2](b) and its inset. The DFT- 
calculated indirect band gap of monolayer ReSe 2 is about 1.22 eV and the direct gap at 
the P point is slightly larger, about 1.25 eV. These are also consistent with previous DFT 
calculationsl^i^i^. 

However, it is well known that DFT usually underestimates the band gap of semiconduc¬ 
tors. Therefore, we have performed the GW calculation to obtain reliable quasiparticle band 
gaps of monolayer ReS 2 and ReSe 2 . Similar to those found in other monolayer 2D semi¬ 
conductors, such as monolayer hexagonal TMDGs (MoS 2 )^^— and phosphorene^, signihcant 
self-energy enhancements are observed in ReX 2 ; at the ’’single-shot” GqWo level, the quasi¬ 
particle band gap of monolayer ReS 2 is increased to be 2.38 eV and that of monolayer ReSe 2 
is increased to 2.09 eV, as listed in Table I. These enhanced many-electron effects are from 
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the depressed screening and reduced dimensionality of suspended 2D semiconductors^‘^>^>^. 

Interestingly, the self-energy not only enlarges the band gap value but also may modify 
the band topology of ReX 2 . This is particularly significant for monolayer ReSe 2 , which 
is an indirect-gap semiconductor at the DFT level, as shown in Fig. [2]^b). However, e-e 
interactions lower the valence band energy at the A point by about 50 meV more, compared 
to that of the F point. Therefore, the VBM is shifted to the F point, making monolayer 
ReSe 2 a direct-gap semiconductor, as listed in the blue-colored VBM in table II. However, it 
has to be pointed out that the indirect to direct transition may not be very sharp because the 
energy difference at the F/ A point are very small and external perturbations may change 
the conclusion. On the other hand, it shall be true that the many-electron self-energy will 
increase the valence band energy at the A point and, at least, make the top of valence band 
almost degenerated. This will affect photoluminescence (PL) experiments. 

These different self-energy corrections are resulted from the nature of the involved elec¬ 
tronic states at the F and A points. For example, a major contribution to the self-energy is 
from the screened-exchange interaction (T^sex)^, which is strongly affected by the spatial 
localization of electronic states. This can be seen from its static form^^: 

EsE.Y(r,r') = -E”i“.^„,(r)C»(r')W'(r,r'), (1) 

where the subscripts n is the band index and k is the sampling points in the reciprocal 
space, (pnk is fhe electronic wave function, and W is the screened Coulomb interaction. 
Obviously, a smaller distance between r and r' enhances the overlap between wave functions 
and screened Coulomb interactions, giving rise to a larger self-energy. This non-uniform 
self-energy correction is similar to what have been observed in hexagonal TMDCs, i.e., the 
different self-energy corrections at the K point and the zone center of monolayer MoS^^. 
Following this idea, we have checked the projected density of states (PDOS) of electronic 
states at the P and A points of ReSe 2 , as listed in Table UTl The A„ valence state has more 
localized d-electron component (86%) than that the P^ valence state (78%). As a result, 
quasiparticle energy of the A.^ state is enhanced more, meaning its energy level is lowered 
more and resulting in a switch of the VBM. 

This direct quasiparticle band gap is not conflicted with experimental measurements, in 
which rather weak PL peak intensity is observed in monolayer ReSe 2 and the PL peak of 
ReSe 2 increases monotonically with increasing the layer number—. This measurement seems 
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to hint an indirect-gap of ReSe 2 . However, different from our calculated suspended and 
neutral cases, these measured samples are on substrates and are inevitably doped. Thus 
their self-energy corrections will be substantially reduced under these conditions, as shown 
in WSe 2 and MoS 2 ^>^, meaning their measured samples may have a slightly indirect band 
gap. It has to be pointed out that this small energy difference is within the intrinsic error 
bars of DFT and GW methods. Thus more accurate experiments of suspended samples are 
important for conclusive results. 

Finally, previous works have shown that the self-consistent GW (sc-GW) scheme be¬ 
yond single-shot calculations may be necessary for 2D semiconductors^>^>^. Therefore, we 
perform one self-consistent update to the Green’s function G; the quasiparticle band gap 
is further increased to 2.82 eV for monolayer ReS 2 and 2.45 eV for monolayer ReSe 2 , as 
listed in Table HI Similar to monolayer M 0 S 2 and black phosphorus, we hnd that more self- 
consistent steps only slightly change the band gap (~ O.leR) and we stop at the sc-GiWq 
level. Finally, the SOG is signihcant in ReX 2 structures (around 150 meV of lowering the 
band gap)^. Gonsidering all these factors, we compute the quasiparticle band gaps and 
summarize them in the last column of Table [B In the following, all our discussions and the 
BSE calculations are based on the hnalized sc-GiWq results with SOG included. 


IV. EXCITON EFFECTS ON OPTICAL ABSORPTION SPECTRA 

Because optical spectra of ReX 2 are anisotropic, we first focus on those with the incident 
light polarized along the F-K2 (x) directions because the corresponding optical absorption 
spectra can exhibit most features. Those of other polarizations will be discussed in Section 

V. First, the single-particle optical absorption spectra for monolayer ReS 2 and ReSe 2 are 
presented in Figs.[3](a) and (c) (blue dashed lines). As expected, the main optical absorption 
shoulder starts from the quasiparticle band gap, which are around 2.7 eV and 2.3 eV, for 
ReS 2 and ReSe 2 , respectively. 

Since e-h interactions are known to be crucial for obtaining reliable optical spectra in 2D 
semiconductors, we have solved the BSE to include them and present the optical absorption 
spectra in Fig. |3]by red solid lines. Similar to other 2D semiconductors, all main optical 
features are dominated by excitonic states in monolayer ReX 2 . For example, after including 
e-h interactions, the hrst absorption peak of ReS 2 is located at 1.63 eV as shown in Fig.[3](a), 
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which is a strongly bound exciton with a 1.07-eV e-h binding energy. Monolayer ReSe 2 
exhibits similarly enhanced excitonic effects as shown in Fig. |3]|^c) with the absorption peak 
at around 1.42 eV and exciton binding energy of 0.87 eV. 

It has to be pointed out that the lowest-energy peak (1.42 eV) in the optical absorption 
spectrum of monolayer ReSe 2 is actually composed of two bright excitons {Ey and F^a), as 
elucidated in the insert of Fig. [31(c). The energy spacing of these two excitons is only about 
96 meV. These two nearly degenerated excitons are resulted from the indirect-gap nature 
of ReSe 2 . To see this, we have plotted the distributions of these exciton amplitudes in the 
reciprocal space (see details in Ref. 0). As shown in Fig. 01 the exciton Ey is mainly from 
the transitions around the zone center (the F point), while the exciton F^a is mainly from 
those around the A points. Since the energies of maxima of the valence band at the F and A 
points are very close to each other, the formed excitons obviously have very similar energies 
as well. 

Recently, several PL measurements of monolayer ReS 2 and ReSe 2 have been performed. 
The measured broad PL peak position is about 1.50 to 1.65 eV at room temperature for 
ReS 2 ^‘^ and 1.47 eV for ReSe 2 ^. These experimental observations are in excellent agree¬ 
ment with our first-principles calculations, which we calculate as 1.6 eV and 1.4 eV for ReS 2 
and ReSe 2 , respectively. As shown in many previous works^^^— , the quasiparticle self-energy 
and e-h interactions partially cancel each other, making the final “optical” gap not much 
larger than the DFT-calculated band gap. 

Finally, we have to address that Coulomb truncation is crucial for all above calculations^IiSS. 
Without Coulomb truncation, both self-energy corrections and excitonic effects will be re¬ 
duced because the actual simulation object is periodic structures with a large interlayer dis¬ 
tance. For instance, for monolayer black phosphorus^, the GW calculation with Coulomb 
truncation gives the band gap as 2 eV while that without truncation gives the value of 1.7 
eV. 

V. ANISOTROPIC EXCITON 

Unlike hexagonal TMDCs, monolayer ReS 2 and ReSe 2 are the distorted 1-T diamond- 
chain with the triclinic symmetry. As a result, the optical absorption may be also anisotropic. 
This has been previously examined at the DFT level^, where it is determined that the ab- 
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sorption spectra are anisotropic with respect to the polarization direction of the incident 
light. However, due to the incorrect band gaps and the short of e-h interactions, those 
DFT-calculations cannot be directly compared with experiments for providing reliable un¬ 
derstanding. This motivates us to study the anisotropic optical responses at the level of 
including quasiparticle energy and excitonic effects. 

Let us hrst focus on the case of monolayer ReS 2 . Its symmetry-induced anisotropic 
optical response can be seen even at the single-particle optical absorption level. As shown 
in Figs. |3l^a) and (b), the single-particle optical absorption spectrum begins at near 2.7 eV 
for the incident light polarized along the F-K2 (x) direction, while it becomes significant at 
2.9 eV for the light polarized along the F-M3 (y) direction. Previous studies have discussed 
this optical anisotropy^*^ and it is attributed to the transitions from non-bonding Re bd 
t 2 g to bd t 2 g and to anti-bonding chalcogen p and a states^. In Table HIl we have listed 
the PDOS of band-edge states of monolayer Res 2 and ReSe 2 . It is clear that both valence 
and conduction band edges are mainly made of Re 5d states and chalcogen p states. For 
example, the F.„ state of ReS 2 has 11 percents of p state of S atoms and 80 percents of d 
state of Re atoms. This is consistent with previous studies^. 

The inclusion of e-h interactions does not essentially change the anisotropic optical re¬ 
sponse. As shown by red solid lines in Figs. [3]^a) and (b), excitonic effects only lower the 
absorption edge but keep the anisotropy. For example, the lowest-energy excitonic peak of 
monolayer ReS 2 is at 1.63 eV for the incident light polarized along the F-K2 (x) direction 
and the absorption edge is located at 1.9 eV for the incident light polarized along the F-M3 
(y) direction. Thus monolayer ReS 2 may work as an optical polarizer covering the frequency 
range from 1.6 eV to 1.9 eV, which is very interesting for near-infrared applications. 

In order to have a complete picture of anisotropies of these important bright excitons, we 
present their optical oscillator strength with respect to the polarization angle of the incident 
light in Fig. |5l As expected, the optical oscillator strength of the prominent exciton of 
monolayer ReS 2 exhibits a £gure-8 shape (Fig.[5](a)), which represents the spatial anisotropic 
optical response. For different directions, the optical oscillator strength can differ by two 
orders of magnitude. For instance, the intensity along F-K2 (x) direction is about 800 
times stronger than that along F-M3 (y) direction. This anisotropic response can easily be 
observed by both optical absorption and PL experiments. 

Monolayer ReSe 2 exhibits a similar anisotropic optical response, with some subtle differ- 


ences. As shown in Figs. [3]^c) and (d), the optical absorption spectra of monolayer ReSe 2 are 
different with respect to the polarization direction of incident light. However, the lowest- 
energy absorption peak (aronnd 1.4 eV) along the r-K2 (x) direction does not completely 
disappear for the incident light polarized along the F-MS (y) direction. This is different 
from the case of ReS 2 . This complication arises from the fact that the lowest-energy peak of 
ReSe 2 actually consists of two excitons {E-p and E\), instead of one in ReS 2 , as shown in the 
insert of Figs. |3l^c) and (d). Ep behaves similar to the lowest exciton in ReS 2 , which is bright 
when the incident light is polarized along the F-K2 (x) direction, and dark along the F-M3 
(y) direction. E\, on the other hand, is not completely dark along the F-M3 (y) direction. 
Thus the lowest-energy peak does not have an extremely high polariztion anisotropy. 

This can be better understood by plotting the oscillator strength of the excitons Ep and 
E\ seperately, which is shown in Figs. EJ^b) and (c). Interestingly, both excitons, Ep and 
F'a, exhibit signihcant spatial anisotropy. However, the preferred direction of the exciton 
E\ is rotated by about 30 degrees. This is not surprising since the main contribution of the 
exciton E^ is from those states at the A points, which are not high-symmetry points. Finally, 
because the prominent absorption peak at 1.4 eV is the combination of these two excitons, 
it is obvious that the overall optical absorption does not show a complete anisotropy as 
shown in Fig. ElJ^d) with a 50-meV smearing. It has to be pointed out that PL experiments 
may exhibit stronger anisotropic optical spectra than the absorption measurements. At 
room temperature, the 96-meV gap between the excitons Ep and i?A will make most excited 
carriers collected by the lower-energy Ep exciton, resulting in a simpler anisotropic PL 
intensity similar to the exciton shown in Fig. El^b). 


VI. SUMMARY 

In summary, hrst-principles GW-BSE simulations have been performed to study quasi¬ 
particle band gaps, excitonic effects, and optical spectra of distorted 1-T diamond-chain 
monolayer ReS 2 and ReSe 2 . Monolayer ReS 2 is direct band gap semiconductor, and the 
quasiparticle band gap is 2.7 eV. Monolayer ReSe 2 changes from the DFT-calculated 1.22- 
eV indirect band gap to be a 2.3-eV direct one after GW corrections. Huge excitonic effects 
dominate optical spectra of suspended monolayer ReX 2 , with around 1-eV e-h binding ener¬ 
gies. The calculated absorption peaks are in excellent agreement with recent measurements. 
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Additionally, both the monolayer ReS 2 and ReSe 2 shows anisotropic optical responses due 
to their reduced triclinic symmetry lattices. In particular, the partially anisotropic optical 
absorption spectra of monolayer ReSe 2 are attributed to the two adjacent excitons, result¬ 
ing from two nearby maxima at the valence band edge. Our calculations provide necessary 
information to clarify experimental measurements and our predictions are useful for broad 
near-infrared applications of distorted 1-T phase 2D semiconductors. 
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TABLE I. DFT and GW band gaps (in unit of eV) for monolayer ReS 2 and ReSe 2 . 



DFT 

GoWo 

sc-GiWo 

sc-GiWo(SOC) 

ReS2 

1.36 

2.38 

2.82 

2.69 


1.43^ 




ReSe2 

1.22 

2.09 

2.45 

2.29 


1.24'^ 





^ Reference IS- 
Reference [l9 . 


TABLE 11. The PDOS of the highest valence and lowest conduction band wave functions for 
monolayer ReS 2 and ReSe 2 at the /c-points labelled in Fig. [2j The percent contributions from each 
atomic orbital for these wave functions are listed. 




Energy (eV) 

DFT GoWo 

s 

Re 

P 

d 

s 

X 

p 

ReS2 

r. 

0.00 

0.00 

0 

5 

80 

3 

11 


Tc 

1.36 

2.38 

1 

2 

69 

1 

27 

ReSe2 

r. 

-0.02 

0.02 

0 

6 

78 

0 

16 


Tc 

1.22 

2.11 

0 

1 

59 

1 

39 


A.V 

0.00 

0.00 

0 

4 

86 

0 

10 


Ac 

1.29 

2.18 

0 

1 

63 

1 

35 
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FIG. 1. (a) Top view of atomic structure of distorted 1-T diamond-chain monolayer ReX 2 . The 
clustering of Re atoms forms the zigzag chains along lattice vector a direction, as shown in the 
red-line box. (b) The side view of monolayer ReX 2 . (c) The first Brillouin zone of monolayer ReX 2 . 



FIG. 2. The DFT-calculated band structures of monolayer ReS 2 (a), and ReSe 2 (b). The inset in 
(b) is the zoom of the band structure inside the red rectangle. The SOG is not included in these 
plots. The top of the valence band is set to be zero. 
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FIG. 3. Optical absorption spectra of monolayer ReS 2 (left) and ReSe 2 (right) for the incident 
light polarized along the r-K 2 direction ((a) and (c)) and the F-Ma direction ((b) and (d)). The 
inset in (c) and (d) zoom in the lowest-energy absorption peaks, which is composed of two adjacent 
excitons {E-p and E\) with an energy difference 96 meV in monolayer ReSe 2 . We employ 50 meV 
Gaussian smearing in these plots. 



FIG. 4. The exciton amplitude distribution of the exciton Ep (a) and that of the E\ (b) (see 
Fig. [3]) in monolayer ReSe 2 in reciprocal space. 
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FIG. 5. Polarization-dependent oscillator strength (arbitrary unit) of excitons in monolayer ReS 2 
(a), the exciton Ey (b) and (c) and their combined oscillator strength (d) in monolayer ReSe 2 
with a 50-meV smearing to mimic the absorption spectrum. 
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Appendix: Band structnre of monolayer ReX 2 with SOC included 


Since Re atoms are heavy atomic species, we thus take into accounts of SOC in Fig. |6l 
Compared with the band structures in Fig. [2J the SOC slightly reduces the band gap of 
monolayer ReX 2 , which is also clearly presented in Table I. 
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FIG. 6. The DFT-calculated band structures of monolayer ReS 2 (a), and ReSe 2 (b) with SOC 
included. The energy of the top of the valence band is set to be zero. 
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